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Abstract 

A new mechanism for an increased event production in the large Q 2 , x region at HERA 
is suggested. This mechanism results from a new type of instanton- induced quark-quark 
interaction, which is related to non-zero quark modes in the instanton field. We estimate 
the contribution of this interaction to the valence quark structure function i 7 ^ using the gas 
instanton model approximation for the QCD vacuum. It is shown that this interaction can 
give a large contribution to F%, especially at large values of Q 2 . The strong dependence of 
this contribution on the quark masses and Q 2 is discussed, and it is concluded that a large 
charm quark contribution to events at high Q 2 can result from charm pair creation in the 
instanton field. We also obtain a sizeable contribution to the valence quark longitudinal 
structure function F£ at high Q 2 . Experimental observables which can check this instanton 
mechanism in the high Q 2 region at HERA are suggested. 



1 Introduction 



Recently[|l|] the HI and ZEUS collaborations at the electron proton collider HERA announced 
an intriguing excess of events above Standard Model expectation in the deep inelastic scattering 
(DIS) region of large Q 2 (> lOOOOGel/ 2 ) and large x. When confirmed by future dataQ, this excess 
is likely to be the most exciting result in particle physics of the last years. Many interesting 
explanations of the this excess, mostly exploring physics scenarios beyond Standard Model, have 
been suggested (see for example ||). 

Within the Standard Model the largest uncertainty on the cross section prediction results 
from the assumed parton densities in the region of large x and large Q 2 . These are based on 
two main ingredients. Firstly, the parton densities are parametrized with a relatively small 
amount of parameters at some low Q 2 scale, using a fit based on pQCD evolution equations 
of experimental data on lepton-parton and parton-parton cross sections at smaller values of 
Q 2 . Secondly, to obtain a prediction for the parton densities at Q 2 > lOOOOGeV^ 2 , leading twist 
(DGLAP) evolution equations || are used. 

A weakness of the first point was pointed out in 0] and ||. In || it was shown that, not 
unexpectedly, the perturbative QCD evolution generates a feed down mechanism, i.e. effects 
at very large x and small Q 2 values propagate to smaller x at larger Q 2 values. Hence e.g. an 
unexpected hard component to the partonic distributions at x > 0.75 and Q 2 ~ 20GeV 2 has 
a visible effect in the region where HERA reports an excess. In || the possible importance of 
an intrinsic charm component in the proton wave function [f|] to explain the excess was pointed 
out. 

However, from our point of view, an important aspect of QCD has not been taken into 
account so far, namely the role of the complicated structure of the QCD vacuum related to the 
presence of strong nonperturbative fluctuations of gluon fields, called instantons ||. In a 
discussion was presented on a possible instanton contribution to the high Q 2 region at HERA. 
It was found that instanton induced interactions can enhance the expected event rate at high 
Q 2 , resulting from multiple emission of gluons at the instanton vertex |10| . 

Instantons describe the tunnelling between different gauge rotated classical vacua in QCD 
field theory, and reflect the nonabelian character of the theory of strong interactions 0. In QCD 
instantons play an important role in chiral symmetry breaking and in the origin of the masses 
of constituent quarks and hadrons. In recent reviews [1^], [13| it was shown that instantons can 



reproduce fundamental quantities of the QCD vacuum, such as the values of the different quark 
and gluon condensates, and can also give a good description of the hadron spectrum. 

Particularly interesting are the interactions which are connected with the so called quark- 



quark t'Hooft |14| interactions, which are induced by instantons. Four important features of 
these interactions should be stressed. 

One feature is the flavor dependence of the interaction, which leads to large contributions in 
flavour singlet channels. For example, the quark-ant iquark repulsion in the t'Hooft interaction 
for SU(3) f singlet qq states leads to a large enhancement of the rf mass and offers a solution to 



the Ua{X) problem [14]. On the other hand, the t'Hooft interaction for quark-quark attraction 
in a flavor singlet state leads to the formation of an isoscalar diquark state inside the nucleon 
| jl5[| . This has direct relevance to DIS: it offers a natural way to explain the experimental data 
on the ratio of proton to neutron structure functions F£(x, Q 2 )/ F^x, Q 2 ) at large x > 0.1, when 

1 with the additional data of the 1997 run the reported excess is somewhat reduced to a two sigma significance 
effect [| 

2 For an introduction to the instanton physics see e.g. [ fTl) . 
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assuming the presence of a rather large isoscalar diquark component in the nucleon, which does 
not disappear for increasing Q 2 WE . 



A second feature of the t'Hooft interaction is it's strong helicity dependence: contrary to the 
perturbative quark-gluon vertex, which conserves the quark helicity, instanton-induced quark- 
quark vertices exhibit a spin flip and consequently lead to quark helicity violation effects. In 
[|17| , fig] ) it was shown that this contribution may explain the observed large violation of the 
Ellis- Jaffe sum rule for the first moment of spin-dependent structure function gi(x,Q 2 ) ||19|| . 

A further property is the strong dependence of the t'Hooft interaction on the mass of the 
quarks involved. For the so called zero quark modes (zero energy states of quarks in the instanton 
field) the contribution behaves as 1/ m* 2 exp(—2m*p c ) , where m* is the effective quark mass in 
the instanton vacuum and p c ~ l.QGeV' 1 is the average instanton size in vacuum||12[. This 
leads to a large suppression of the contribution of zero quark modes to the quark distribution 
function for heavy flavours. A simple estimation in the framework of the instanton liquid model 
1^] shows that sea charm quarks are suppressed in F 2 by a factor 10 -4 compared to light sea 



quarks, and can therefore be neglected. 

A final feature is related to the strong dependence of the t'Hooft interaction on the virtuality 
of the interacting quark. The t'Hooft interaction results from the existence of zero quark modes 
in instanton fields, hence they are strongly localized in the field, which leads to an exponential 
dependence of their Green function on the quark virtuality. 

One can calculate the instanton contribution to the sea quark distribution function by con- 
sidering the interaction of a quark from a virtual photon with a valence quark in the nucleon 
through an instanton fluctuation (see Fig. 1). In this case all properties of the instanton contri- 
bution are determined by the dependence of the instanton-induced quark-quark cross section on 
the virtuality of the incoming sea quark. Using the model of Landshoff [20] one can show that a 
fast reduction of the strength of this interaction with increasing sea quark virtuality results into 
a leading twist contribution in the sea quark distribution function [|l7j] (see also pll ). 

The exponential dependence of the quark-instanton vertex for zero quark modes on the in- 
coming quark virtuality is expected to lead to a rather soft instanton-induced quark distribution 
as function of Bjorken x. This makes it difficult to explain a substantial excess at large values 
of x in the HERA region from quark zero modes in the instanton field (see ||) alone. 

The main goal of this study is to consider the contribution of so called non-zero quark modes 
in the instanton field (i.e. non-zero energy states of quarks in the instanton field) to the F 2 (x, Q 2 ) 
and Fl(x, Q 2 ) structure functions, with special emphasis to the large Q 2 and x region at HERA. 
In most of the instanton calculations performed so far the contribution of non-zero quark modes 
has been neglected. It has been argued that this contribution is proportional to the product 
m q p c , which is indeed small for light u and d quarks. However, e.g. for strange quarks, this value 
amounts to m s p c « 0.25. Evidently for the heavy quarks c, b and t, the contribution of non-zero 
quark modes in the instanton field will dominate over the contribution from zero modes. 

We show that the contribution of instanton interactions to unpolarized structure functions 
for these modes become increasingly important at large Q 2 . Hence they can provide a QCD 
mechanism which leads to an any excess seen above calculations using parton distributions 
which result from 'standard' perturbative QCD evolution. For a quantitative prediction of these 
contributions we will use gas instanton model approximation for the QCD vacuum. In section 2 
we calculate the instanton induced contribution to the quark structure functions, using the Green 
functions formalism. In section 3 we study the structure of the interaction using the effective 
Lagrangian formalism. The calculation of the non-zero modes contribution in framework of the 



liquid instanton model JT2L |T3| will be the subject of a future study. 
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2 Instanton induced contribution to the quark structure 
functions 



To estimate the instanton induced contribution from non-zero quark mode states to the 
nucleon structure functions F 2 (x, Q 2 ), F L (x,Q 2 ) we consider the contribution to valence quark 
structure function which is given by the diagrams in Fig.l. 





a] 



b) 



Figure 1: The contribution from sea quarks created by instantons to the valence quark structure 
functions, a) the interaction of the sea quark with the valence quark and b) the interaction of the 
sea antiquark with the valence quark. 



This contribution is determined by the deviation of the quark Green function in the instanton 
field (i.e. the quark propagator in the background colour instanton field) from a free quark Green 
function Sq(x, y). 

S mt (x,y) = S I (x,y)-S (x,y), (1) 

with x,y Euclidean space coordinates, and where Si(x,y) is the quark Green function in an 
instanton field and 

a i \ 1 • ( x ~ V) fo \ 

S °^ y) = - 2^(x-yr (2) 

The function Si(x,y) for a quark with current mass m can be presented as the sum of the zero 
quark modes and non-zero quark modes contribution (see |12]) Q 



Sr(x, y) = S z (x, y) + S nz (x } y) = — + ^ 



im 



A + im 



(3) 



3 We will present most formulae in Euclidean space. Only in the final step of the calculation the analytical 
continuation to the Minkowskii space will be performed. 
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where 

*oiA x )*ojAy) = g^MyK^ivyi-h ® {Ur-T+u) a p, (4) 

<f>{x) = (5) 

VT v ^2( X 2 +p 2)3/2 

and 

q / x 1 1 IQ , W1 , p 2 Ur~ ■ xr + ■ yU+ 



p 2 Ut XT+--/T -(x-y)T + -yU + 
D (x, yj^-^i y~ — 7+ + 

Ut~ ■ xt + ■ (x — y)r~ ■ 7t + • yU + 



p 2 + y 2 



7 _)) + 0(m), (6) 



where i,j (a, (3) are the spinor (colour) indices, j± = (l±7 5 )/2, = (r, =f) is the colour matrix, 
U is the orientation matrix of instanton in the colour space, p is the instanton size, and 

D ^ V) = 4*2(x- yf (7) 

is the propagator of a scalar quark. 

The contribution of zero quark modes to the polarized structure function g\{x) and the 
unpolarized structure function F 2 has been estimated in 0] and 0. It was shown that this 
contribution is rather large, specially at low x, can explain the violation of the Ellis- Jaffe sum 
rule and gives a contribution to the high Q 2 region at HERA, on top of the perturbative QCD 
contribution. 

However, as discussed in section 1, in the large x and Q 2 region one can expect also large 
contributions from non-zero quark modes, which are not very strongly localized in the instan- 
ton field and therefore their contribution is expected to have a different x and Q 2 dependence 
compared to the one from zero quark modes. 

We will estimate this contribution from non-zero quark modes for the case of two flavor quark 
Nf = 2, where one quark will be heavy. In this case the dominated contribution is related to 
the one shown in the diagrams of Fig.l, where the sea quark is in a non-zero mode state and the 
valence quark is in a zero mode state. 

The contribution to the valence quark part of the structure functions F% (a;, Q 2 ) and F[(x, Q 2 ) 
is given by 

F"(x,Q 2 ) = (- gtlv + 6x^)xW,„(q,p), 

p-q 

F[(x,Q 2 ) = Ax 2 ^W, u (q,P), (8) 

where p is the momentum of the valence quark, q is the virtual photon momentum, x = Q 2 /2p-q 
and 

W,u(p,q) = ^f[ J ^(27r)%k 2 -m 2 )5(q + p-h-k 2 -k 3 )T,(p,q)T:(p,q), (9) 



where T M (p, q) is the matrix element of the virtual photon- valence quark scattering. 
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The instanton-induced non-zero quark mode contribution in leading order in the mass ex- 
pansion, for sea quarks with mass m q , is 

T>, q) = -ie q fj" J dSe^Kh)h{S**{y> k 2 ) llM S nz (h, y) + 

^(y.^T^^i^))*!*^)^)^^^^)^), (10) 

where s(ki) and u(p) are Dirac spinors for sea quarks and valence quark respectively, and p cu t is 
a parameter which cuts the contribution from large instantons. The density of the instantons is 

d( P ) = ^ e -**+"/<* [-^)e~^{ppt flKp), (11) 

where 

b = b + ^(b 1 -l2b ), (12) 



a s (p) 



An 



&i lnln(^ 



and in the MS-scheme we have Ci = 0.466, C 2 = 1.54 and C3 = 0.153, where \i is the renormal- 
ization scale, and 

60 = 11-^ b 1 = 102 -jN f . (13) 

For numerical estimations below the value A = 230MeV was taken. 

The Fourier transform of the Green function for non-zero mode states for quarks is 



S nz (y,h) = j d 4 xe lk2X S nz (y,x) 
and antiquarks is S nz (ki, y) for the case k\ — > 0, fcf — ► has been calculated in pi 



(14) 



k\S nz {k u y) 



2y 2 



^2 • 2/ 



nj 



Vy 2 + p 2 



Jkiy 



1 I P 2 n U T--kiT + -yU + (1 c _ iklV , 



2y* 



h ■ y 



(15) 



The Fourier transformed Green function for zero quark modes in an instanton field equals 



p 2 7T 2 



hS z {h,p)p = *L_( V7 „ 7 _) ® (Ur^U + )F(p 2 )F(kl), 



where 



F(k 2 



-> 1 at k 2 -> 0. 



(16) 
(17) 



The result of the calculation of fllP]), using (|15]), (|16|) at kf — > is 



-e 9 7r 



p 



ufaXj^j-) ® (f/ V r+f/+) M (p) 



m,, 



;is) 



where k = q — k\ — k 2 and 



, , , Ut~ ■ k 2 r + ■ kU + 



F T (k, k lt k 2 ) = h(k 2 ) + ft— I 2 (k,k 1 ,k 2 ) + 

k 2 ■ k 

Ur~ ■ k 2 r+ ■ hU + r n . Ut- -kr+ ■k 1 U+ _ n , , 
-h{k, ki, k 2 ) H I±(k, h, k 2 ) + h 



with 



k 2 ■ k\ " " k\ ■ k 



2 = J_ r° dzZ 3 (z-y/z*-p?W) 

p 2 k A Jo z 2 — p 2 k 2 



I 2 (k,k u k 2 ) = p 



I s ik : ki,k 2 ) = p 2 k 2 -ki [ dp 

Jo 



~{k + k 2 ) 2 V^k 2 
i KoU-ik + h + pktfp) - K {J-(k + /3k 2 ) 2 p) 



2(k ■ k x + phi ■ k 2 ) 



vw^mr (z) - (19) 

To get the final result we have to integrate over the colour orientation dU which is given by 

dUUijU^ = —5 jk 5 u 

I i 

dUUijU M U mn U op = j-j^SjkSudnodmp + ^jy 2 ^ji^Jjki^Jlii^ )no(A )mp- (20) 

In this paper only the colour singlet exchange between quarks induced by instantons is taken 
into account. The contribution of colour octet exchange is suppressed by a factor 1/N C and will 
be studied in a forthcoming paper. For the colour singlet case the final result for the non-zero 
quark modes contribution to the valence quark structure functions F 2 and Fl, after averaging 
over the initial quark polarization and colour, equals 

where dPS 3 is the phase space integration, t 2 = (p — fc 3 ) 2 , and 

Fj{k, h, k 2 ) = h{k 2 ) + I 2 (k, h, k 2 ) + I 3 (k, h, k 2 ) + h(k, h, k 2 ) + k x ^ k 2 , (22) 

The integration is performed in the center of mass frame of the sea quarks. In this frame the 
momenta of particles are (see for example |[22|| ) 



q = (Q ,0,0,E q ) 

t = (t Q ,0,0,-E q ) 

p = E p (l, sin/3, 0, cos/5) 

k\ j2 = (Ek, TQk sin 9 cos ^fq k sin 6 sin $, cos#), (23) 
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where 



Ek 
Qo 



1/2, q k = yj8- 4m 2 /2 
s-Q 2 - t 2 )/4E k , t = 2E k - Qo, 



„ v ( S + Q2 + i 2)/ 45 _ f2j 5 G =(g + p) 2 , 
£ P = (S G + Q 2 + t 2 )/4£ fc , s = (t + g) 2 , 



cos /? 



- Q 2 + (S G + Q 2 + t 2 )Q /2E k )/(2E p E qJ 



For the phase space integration we have 

x r s c 



dPS d 



ds 



dt 2 



dt\ 



d<t> 



1287T 3 Q 2 75 A) ' h\ min s-t 2 + Q 2 27r' 
where t\ = (q - h) 2 , S = (2m q ) 2 , S G = Q 2 (l - x)/x, t 2 max = -(S G - s)/(l - x), and 



(24) 



(25) 



jiiiax.min 



s + Q 2 -t 2 



1 ± 



\ 



(i-^)(i+ 



iQ'V 



's + Q 2 -t 



2\2- 



(26) 



In the calculation we take for the renormalization scale fi = l/p cut . For most calculations 
there is a very strong dependence on the cutoff scale p cu t of the instanton density on instanton 
induced contributions to physical quantities. From our point of view a natural way to derive 
an estimate for this parameter is by considering the confinement of the valence quarks inside 
the proton. For a zero mode valence quark in the instanton field we thus expect a cutoff factor 
exp(— 2\p\p), where \p\ ~ 300MeV represents the average virtuality of a valence quark in the 
proton. This leads to a value of p cut ~ 1.6GeV _1 , which happens to coincide with the average 
size of the instantons in QCD vacuum (see [0). 

In Fig. 2 we present the result of the calculation of the instanton induced contribution from 
the non-zero charm quark modes to the structure functions F%(x, Q 2 ) and F£(x, Q 2 ) at x — 0.5. 
This is approximately the x region where the HERA experiments have reported an excess at 
high Q 2 . In the gas instanton model and for case Nf = 2, as studied here, the strange quark 
contribution to the F£ and F£ structure functions from non-zero quark modes is expected to 
be about a factor 10~ 2 suppressed compared to the one of the charm quark, due to the smaller 
current quark mass. However, in the framework of the perhaps more realistic liquid instanton 
model for the QCD vacuum [|l!J, [|3| an additional enhancement of strange quark production 
is expected, due to chiral symmetry breaking which leads to an increase of the effective quark 
mass. 

In Fig. 3 the result of the calculation of instanton contribution to the F2(x,Q 2 ) structure 
function is compared with the perturbative QCD contribution to the same structure function, 
due to gluon emission 



{x,C? 



2tt 



-p qq (x)H 



Q 2 



where 



4 1 + x 2 
31-i 



(27) 



(28) 
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Q 2 GeV 2 



Figure 2: The instanton-induced charm quark pair contribution to the valence quark structure 
functions F^x^Q 2 ) at x = 0.5 (full line) and Fl(x,Q 2 ) (dashed line), at x = 0.5, for a charm 
quark mass of m c = 1.5GeV. 

The instanton contribution to the parton distributions appears to be as large as as the perturba- 
tive evolution correction from gluons. This suggests that this effect should be taken into account 
in the Q 2 analysis of the proton F$(x, Q 2 ) structure function at a large x and Q 2 . 

The result for the instanton induced non-zero quark modes contribution is proportional to 
the square of the mass of the quark which is in the non-zero mode state. In our calculation 
we have taken only into account terms of order m 2 , but we did not explicitly take the sea 
quark mass in the Green function into account. At small Q 2 ~ 20GeV 2 this contribution is 
rather small due to phase space restriction. But for large Q 2 the contribution to both F£(x, Q 2 ) 
and Fl(x,Q 2 ) becomes very large due to the strong Q 2 dependence of this contribution. This 
strong Q 2 dependence of the instanton contribution originates from the dimensional effective 
coupling of the instanton-induced interaction which is given by some power of the instanton 
size p. Therefore, similar to the famous Fermi effective weak interaction, the interaction through 
instantons should in principle lead to a violation of the unitarity bound for the DIS cross section. 
However, additional perturbative and nonperturbative QCD corrections to the bare instanton 
vertex, as discussed in in |23], , are expected to cancel the anomalous growth of the instanton 
contribution to DIS structure functions at very large values of Q 2 . 

It is interesting to note that the non-zero quark modes also lead to a sizable contribution 
to the Fl(x,Q 2 ) structure function at large Q 2 . This results from the structure of the Green 
function for quarks in the instanton field which includes a part that comes from to the Green 
function of scalar quark (see (|6]), (|7|)). 

Therefore the charm quark contribution induced by instantons to structure functions is large. 

4 Threshold effect due to the quark masses are taken into account in the phase space integration. 
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Figure 3: The ratio of the instanton induced charm quark pair contribution to the valence quark 
structure function F^x^Q 2 ), to the perturbative QCD contribution, at x = 0.5. 

We also expect some contribution from b quarks due to the same mechanism. Experimentally 
this effect could be verified by studying the evolution of the charm component of the structure 
function. A considerable increase with increasing Q 2 is predicted. 

While this study was obviously stimulated by the excess of the high Q 2 events seen at HERA, 
it should be noted that independently of it's confirmation in future, this type of contribution 
can play an important role in QCD analyses of very high Q 2 data, and should therefore be 
considered. In a future paper we plan to make more complete study of its implications for QCD 
studies. 



3 Effective quark-quark interaction induced by non-zero 
quark modes in instanton field 



In the previous section a striking Q 2 dependence of the non-zero quark modes contribution 
to the structure functions has been found. This Q 2 dependence results from the structure of 
quark-quark interaction vertex induced by non-zero quark modes in the instanton field. Let 
us consider an instanton induced two quark interaction for the case when one of the incoming 
quarks is off-shell, a situation which is directly related to the sea quark contribution to the 
structure function. One part of the interaction is related to the t'Hooft interaction induced by 
quark zero mode states. The effective vertex can obtained by amputating of the four-point zero 
modes quark Green function in instanton field 



I 




u(k 2 )k 2 S z (k 2 , ki)kiu(ki)d(k 3 )k 3 S z (k 3 , p)pd{p) . 



(29) 
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After substitution the Fourier transformed zero modes (|T^) and performing the integration over 
the instanton orientation in colour space dU, one obtains the following effective Lagrangian for 
the t'Hooft interaction 



2^2 



£.= [ ^^(^) 2 (%(^2)^(fci)4(A:3)rfL(p)(l+^(l-^>^J«+(i? M L)))F(|A; 1 |p), 

(30) 

where F(\k\\p) m exp{—\k\\p) for p — > 0. The density of the instantons is proportional to the 
product of the quark masses ([□]) and therefore, in spite of the quark masses in the denominator 
of Eq. (0), the final result is finite in the limit m u , — > 0. 

The effective Lagrangian connected with the non-zero modes contribution to leading order 
in m q , and for a small instanton size p — > can be obtained in a similar way 

l[ f> s(k 2 )k 2 {S nz (k 2 , h) - S (k 2 , h)} k\s{ki ) d(k 3 )k 3 S z (k 3 ,p)pd(p). (31) 



P 5 



The final result is 



32 k 2 ■ k 1 
where t 2 = (k\ — k 3 ) 2 and the form factors 



dpd(p) ^ 2 p\ 2 {„ , 3^ aw _^ , 
m d p 5 

d Ff(klt 2 ))® sikJhsikMhMp), (32) 



9 t a s t a d k T 2 k1a. d ^ 





Fm,e) = ,.- == - ¥WTw )Jl(z) . (33) 

A very interesting feature of this interaction (|32|) is that it is zero for on-shell massless quarks 
but it is large for the case where one of the quarks is off-shell. This is the origin of the striking 
Q 2 dependence of the contribution of this interaction to the DIS structure functions, as found in 
the previous section. Another feature of the non-zero mode interaction, as compared to the zero 
mode interaction ( |3~TD is the helicity flip. The t'Hooft interaction leads to a double spin flip of 
quarks and therefore this interaction can give a large contribution to the double spin asymmetry 
in the polarized DIS (see |17|]). The contrary is true for non-zero quark mode interactions, which 



do not have a definite helicity and can lead only to a single spin flip; their contribution to the 
double spin asymmetry is zero . Hence these instanton induced interactions do not contribute 
to the spin-dependent structure function g±(x, Q 2 ). The zero quark mode interactions on the 
other hand do contribute to g\(x,Q 2 ), and this contribution is negative 0. Thus the total 
contribution of instantons to gi(x, Q 2 ) is negative, contrary to the perturbative QCD expectation 
which yield a positive value for g±(x, Q 2 ) at large x and Q 2 . Hence the instanton mechanism for 
DIS events excess at HERA can be checked by a measurement of the double spin asymmetry 
for these events at high Q 2 , when a polarized proton beam at HERA will be available. The zero 
quark mode component will contribute with a negative asymmetry, while the non-zero quark 
mode is predicted to reduce the asymmetry, compared to the positive asymmetry predicted from 
perturbative QCD. 
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4 Summary 



In summary, a new contribution to the proton structure at high Q 2 and high x is suggested. 
This contribution is connected with complicated structure of the QCD vacuum, related to the 
existence of strong fluctuations of the vacuum gluon fields: instantons. It is shown that a new 
type of quark-quark interaction in QCD, which is due to non-zero quark mode states in the 
instanton field, gives a large heavy quark contribution to the quark structure functions in the 
large Q 2 , x region. This effect could already be visible at HERA. 

The large Q 2 dependence of this interaction comes from the dependence of the effective 
instanton-induced vertex on the quark virtuality. We also predict anomalous spin properties of 
high Q 2 events induced by this new interaction. Due to the strong dependence of this interac- 
tion on the quark virtuality it should show up also in other lepton-hadron and hadron-hadron 
processes with large transfer momentum. Note that this contribution is missing in standard 
QCD analyses. Independently of a possible excess at large Q 2 and x of the HERA data this 
contribution may have to be be considered when evolving to large scales. Finally we note that 
enhancement of the charm contribution to the structure function at large x can lead also to an 
explicit QCD model for intrinsic charm ]7| in the proton. 
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